home *** CD-ROM | disk | FTP | other *** search
Text File | 1993-09-19 | 468 b | 25 lines | [TEXT/MPAD] |
- -- Solution of a system of equations by Cramers rule.
-
- solve(A,B)[j] = det(cram(A,B,j))/det(A)
- cram(A,B,k)[i,j] = A[i,j] when j≠k, B[i] dim[count(B)]
-
- -- example system of equations
- -- 2x + y + z = 2
- -- x - y + 5z = 4
- -- y - z = 4
-
- A={{ 2, 1, 1}, -- coefficients
- { 1,-1, 5},
- { 0, 1,-1}}
-
- C= { 2, 4, 4} -- constants
-
- solve(A,C):{-4.0,7.0,3.0}
-
- x=solve(A,C)[1]
- y=solve(A,C)[2]
- z=solve(A,C)[3]
-
- 2*x + y + z:2.0
- x - y + 5*z:4.0
- y - z:4.0